Evidence for the ground-state resonance of ^^O 
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Evidence for the ground state of the neutron- unbound nucleus ^^O was observed for the first 
time in the single proton- knockout reaction from a 82 MeV/u ^^F beam. Neutrons were measured 
in coincidence with fragments. ^^O was determined to be unbound by ISO^iso keV from the 
observation of low-energy neutrons. This result agrees with recent shell model calculations based 
on microscopic two- and three-nucleon forces. 



A major challenge in nuclear physics remains the de- 
scription of nuclei based on fundamental interactions. 
"Ab-initio" approaches have been developed to calculate 
nuclear properties based on nucleon-nucleon scattering 
data up to A ^ 12 [T]. Recent advances in nuclear the- 
ory made it possible to describe some fundamental prop- 
erties of light nuclei up to oxygen based on two- and 
three-nucleon interactions [2HS]- On the way to heavier 
nuclides it will be critical for these models to describe the 
dramatic change in the location of the neutron dripline 
from oxygen {N — 16) to fluorine [N > 22) which was 
first pointed out by Sakurai et al. [7j. The addition of 
one proton binds at least six additional neutrons. The 
two- neutron separation energy of ^^O serves as an im- 
portant benchmark for these calculations. The majority 
of the current nuclear structure models predict ^^O to be 
bound [HHlSj. Experimentally it has been shown that ^^O 
is bound [Tl] while repeated searches for bound ^^O and 
^^O have been unsuccessful [TSHII], although ^^O had ini- 
tially been reported as being particle stable [22] • Shell- 
model calculations using phenomenological interactions 
do predict to be unbound: SDPF-M [23] by 77 keV 
and USDOSa [24] by 510 keV. A continuum shell model 
calculation predicts ^^O to be unbound by only 21 keV 
[25j . Recently it was shown that three-body forces are 
necessary to describe the binding energies of neutron-rich 
oxygen isotopes based on fundamental nucleon-nucleon 
forces However, no calculations have been published 
which simultaneously predict ^^O to be unbound and '^^F 
to be bound. 

Because no bound states of exist, the search for 
its elusive ground state must be extended to unbound 
states. The unbound ground state of ^^O was measured 
using invariant mass spectroscopy and was found to have 
a decay energy of 7701^q keV [26] • Due to this high 
ground-state energy of it is hkely that ^^O is bound 
with respect to one-neutron emission and unbound with 
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FIG. 1: Schematic decay level scheme of O. Proton- 
knockout reactions from ^®F populate the ground state of ^^O 
which was measured to decay to the ground state of (solid 
arrow and lines) [26]. The knockout reaction from ^^F used 
in the present work populates states in ^^O. The dashed lines 
show the predicted levels calculated by the continuum shell 
model [25^. Possible decay channels from ^^O are shown by 
the dotted arrows. 



respect to two- neutron emission. ^^O is thus also an ex- 
cellent candidate for di-neutron emission. Furthermore 
calculations by Grigorenko et al. predict that the emis- 
sion of a pair of correlated neutrons might be hindered 
so that for very low decay energies lifetimes on the order 
of pico- to nanoseconds could be possible [57] ■ 

We searched for unbound states in ^^O using one- 
proton knockout reactions from ^''F and by measuring 
neutrons in coincidence with ^^O fragments. Figure [l] 
shows a schematic level scheme of the possible decay 
paths for predicted states of ^^O. In this letter we present 
the first evidence for the observation of the unbound 
ground state of ^^O. 
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Particle identification spectrum. The three oxygen 
^^O, ^^O, and ^''O were identified as indicated in 



FIG. 2: 

isotopes 

the figure. The grey-shaded area corresponds to the events 
included in the present analysis for the decay of ^^O and ^^0. 



The experiment was performed at the Coupled Cy- 
clotron Facility of the National Superconducting Cy- 
clotron Laboratory at Michigan State University. A 
primary beam of 140 MeV/u ''^Ca^""'" bombarded a 
1316 mg/cm^ ^Be production target. The desired ^^F 
secondary beam with an energy of 82 MeV/u was sep- 
arated from the other reaction products, and primary 
beam, using the A1900 fragment separator [5H1 with a 
1050 mg/cm^ Al wedge placed at the intermediate fo- 
cal plane. The ^''F component of the secondary beam 
was 7% with the main contaminant being "^"Na. The ^^F 
fragments were identified event-by-event through time-of- 
flight, and the average intensity was 14 per second with a 
momentum spread of 2.5%. The ^''F beam then impinged 
on a 705 mg/cm^ Be target producing the isotope of in- 
terest, ^^O, through one-proton removal reactions. Mea- 
surement of the ^^O — > + 2n decay required detec- 
tion of both the neutrons and a charged particle. Beyond 
the reaction target, the large-gap 4 Tm superconduct- 
ing dipole (Sweeper) magnet [29] was used to bend the 
charged particles 43° after which they passed through 
a set of position and energy sensitive detectors. The 
fragments were identified by time-of-flight and en- 
ergy loss in the charged particle detectors. The Modular 
Neutron Array (MoNA) [30] was placed 6.05 m from the 
reaction target and measured the angle and energy (from 
time-of-flight) of the forward-focused beam- velocity neu- 
trons. The charged particle detectors and MoNA provide 
the relativistic four-momentum vectors for the ^^O nuclei 
and neutrons that were used to calculate the decay energy 
of the two-body (fragment -I- n) or three-body (fragment 
+ n + n) systems. The thicker target, as compared to 
previous experiments, increases the FWHM of the decay 
energy resolution to 200 and 800 keV for decay energies of 
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FIG. 3: Decay energy spectrum of ^^O. The data points 
are from the present one-proton knockout reaction from ^^F, 
while the solid line corresponds to the fit to the data of the 
one-proton knockout reaction from ^®F from reference [26) . 



100 and 800 keV, respectively. Further details about the 
experimental set-up, parameters and analysis procedures 
can be found in references [5S1 [5TH5i] . 

Figure [2] shows the particle identification spectrum for 
the oxygen isotopes. The events included in the further 
analysis are shown by the grey-shaded area. With this 
stringent cut the contamination of ^"^O events was limited 
to about 1%. 

The decay energy spectrum for the two-body (fragment 
+ n) system corresponding to ^^O is shown in Figure [s] 
^^O previously had been measured and exhibited a reso- 
nance with a decay energy of 770t lo keV |26j . The solid 
line is representative of this resonance. It is immediately 
obvious that the data from the present one-proton knock- 
out reaction from ^^F (data points) has a completely dif- 
ferent shape with an additional strong peak at low de- 
cay energies. Because these low-decay energy neutrons 
cannot come from ^^O they have to originate from ^^O. 
This initial evidence for a low energy resonance in ^^O 
is substantiated by the presence of a low-energy peak in 
the reconstructed three-body (fragment -I- n -I- n) system 
corresponding to ^^O as shown in Figure [4](b). It should 
be mentioned that in the present experiment as well as in 
the experiment of reference [34 data were also recorded 
for neutrons in coincidence with ^^O -I- n, and the mea- 
sured decay energy spectra from the two experiments are 
in agreement. 

The three-body decay energy spectrum shown in Fig- 
ure[4](b) was reconstructed from a ^''O fragment and two 
neutron interactions in MoNA. Thus, the spectrum also 
contains cross-talk events from a single neutron scatter- 
ing twice. The broad distribution of counts up to approx- 
imately 3 MeV in the three-body decay energy spectrum 
could be due to the ground state of ^^O where the neu- 
tron scattered twice and/or the sequential decay from 
excited states in ^^O through the ^^O ground state. In 
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FIG. 4: (Color Online) (a) Decay energy spectrum of ^^O (two-body, fragment + n system). The data are the same as in 
Figure Isl The inset indicates the different decay paths. (b,c,d) decay energy spectra of (three-body, fragment -I- n -I- n 
systemJT^ The grey area in (b) represents the simulated cross-talk contribution. The individual contributions to the cross-talk 
from the different simulated states are shown in the inset. Panel (c) shows the cross-talk subtracted spectrum. Panel (d) 
shows the spectrum with causality cuts applied to the data as well as the simulation. The lines are explained in the text. 



order to fit the data, we performed Monte Carlo simula- 
tions which included the geometrical acceptances, energy, 
positions, and timing resolution, tracking of the charged 
particles through the Sweeper magnet, and the reaction 
and decay mechanisms. The interaction of the neutrons 
with MoNA was described using the Geant4 simulation 
package [35] with the addition of the MENATE-R physics 
class |36] . Thus, multiple interaction of a single neutron 
were fully included in the simulations. 

A low-energy resonance was simulated with a Breit- 
Wigner line shape and the energy from the ^^O decay 
into ^'*0 and two neutrons was partitioned between the 
three outgoing particles according to the phase-space 
model of references [371 [3H] ■ The data are not sensitive 
to the detailed parameters of the various possible con- 
tributions to the high energy continuum. As shown in 
Figure [ij the continuum shell model predicts only one 
excited state for ^^O at approximately 2 MeV with the 
next excited state calculated at about 6 MeV Thus, 
a resonance in ^^O at a fixed energy of 2 MeV with a 
width of 200 keV, which was allowed to decay sequen- 
tially via the known unbound ground state of ^^O, was 
included in the simulation. The resonance parameters for 
'^O (£;docay = 770 keV, Fdecay = 172 kcV, L = 2) were 
taken from reference [2S]. A fit to the two-body and 



three-body system was performed where the resonance 
energy and width of the low-energy resonance and the 
relative strengths of the three contributions (low-energy 
and 2 MeV state in ^^O and the ground state of ^^O) 
were free parameters. 

Figure |4] shows the resulting simulated spectra for the 
best fit parameters (solid line) . The two-body and three- 
body systems are shown on panels (a) and (b), respec- 
tively. The low-energy resonance in ^^O at i^decay — 
150 keV and Fdecay = 5 keV is shown by the long-dashed 
red line and the 2 MeV resonance by the short-dashed 
green line. The dot-dashed blue line shows the contribu- 
tion from the direct population of the ^^O ground state. 
These decay paths are also indicated in the inset of Fig- 
ure |4ja) . In the simulation it is possible to distinguish 
real two-neutron detection from cross-talk events where 
a single neutron interacted twice in MoNA. The grey- 
shaded area in Figure [4] (b) shows the contributions of 
the simulated cross-talk events to the total spectrum. In 
the inset the individual contributions to the cross-talk 
from the low- and high-energy ^^O decays as well as the 
contribution from the ^^O are shown. In order to demon- 
strate the positive signal of real 2n events we applied 
two different methods. First we subtracted the simu- 
lated cross-talk events from the data, and the results are 
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shown in Figure ^c) . While the high-energy events are 
essentially reduced to background level, the low-energy 
peak is clearly still present. In a second method we ap- 
plied causality cuts to the data. Following the description 
of references [321 EO] we required a spatial (25 cm) and 
velocity (7 cm/ns) separation of two interactions in the 
data as well as in the simulation. At the expense of a re- 
duction of ~50% in efficiency, the cross-talk contribution 
was reduced by a factor of three. The results are shown 
in figure gd) which again clearly shows the presence of 
the low energy peak of real two-neutron events. 

The best fit to the data included a resonance for the 
^^O ground state of ISOl^gg keV. This value agrees with 
the recent calculations of a low-energy unbound reso- 
nance in ^^O [11I2S]. The fit was insensitive to the width 
of the resonance. The cross section for populating the 
^^O ground state was 1.8±1.0 mb while the cross section 
for populating ^^O was 4.2±2 mb. These values are con- 
sistent with Ip and Ipln removal cross sections from ^''F, 
^^F and ^^F. While the measured values for ^''F agreed 
with the calculated removal cross sections the cross sec- 
tions decreased with increasing neutron number |18j . 

One might speculate that the low-energy neutron orig- 
inates from a state located around 900 keV and then 
decays with a ^100 keV neutron sequentially via the 
ground state of ^^O. The observed low energy peak would 
then be solely due to the low-energy neutron scatter- 
ing twice in MoNA. We simulated such a decay and the 
results are shown in Figure [5] In addition to the to- 
tal fit (solid line), the figure also shows the individual 
contributions from real two-neutron coincidences (long- 
dashed blue line) and from events where either the low- or 
the high-energy neutron scatters twice in MoNA (short- 
dashed red line). Clearly, this sequential decay scenario 
is not supported by the data. 

In conclusion, we present evidence for the observation 
of the ^^O ground state which is unbound by less than 
200 keV and which decays by emitting two low-energy 
neutrons. A future experiment with higher statistics is 
necessary to study the detail of the decay mechanism, 
i.e. explore the possibility of a di-neutron decay. The 
upper limit for the total decay energy leaves open the 
exciting prospect that a di-neutron decay of ^^O might 
have a long half-life (> picoseconds) as first suggested by 
Grigorenko et al. |27| . 
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FIG. 5: (Color online) Decay energy spectrum from two 
neutron interactions in coincidence with ^^O (same as Figure 
|4]^b)). The solid line corresponds to the results of a simu- 
lation assuming a hypothetical 870 keV state in '^''O decay- 
ing sequentially via the ground state of '^^O. The dashed and 
dot-dashed lines are the contributions due to real two-neutron 
coincidences and double-hits from one neutron, respectively. 
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